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Local Spectral Density and Vacuum Energy 
Near a Quantum Graph Vertex 

Stephen A. Fulling 

Abstract. The delta interaction at a vertex generalizes the Robin (general- 
ized Neumann) boundary condition on an interval. Study of a single vertex 
with N infinite leads suffices to determine the localized effects of such a vertex 
on densities of states, etc. For all the standard initial-value problems, such as 
that for the wave equation, the pertinent integral kernel (Green function) on 
the graph can be easily constructed from the corresponding elementary Green 
function on the real line. From the results one obtains the spectral-projection 
kernel, local spectral density, and local energy density. The energy density, 
which refers to an interpretation of the graph as the domain of a quantized 
scalar field, is a coefficient in the asymptotic expansion of the Green func- 
tion for an elliptic problem involving the graph Hamiltonian; that expansion 
contains spectral/geometrical information beyond that in the much-studied 
heat-kernel expansion. 



Introduction 

A topic of perennial and renewed interest in quantum field theory is the en- 
ergy of the "vacuum" — that is, of the ground state of a field subjected to some 
nontrivial external condition [121 [13 El HH IHl |H1- (The prototype is the 
electromagnetic field between two parallel flat conductors.) Although the only 
quantities indisputably open to experiment are the derivatives of the total energy 
with respect to parameters defining the configuration, the energy itself and even 
its localization in space are of theoretical interest, not least because energy den- 
sity (along with associated quantities such as pressure) acts as the source of the 
gravitational field in general relativity |23l I18| . 

From a mathematical point of view, vacuum energy is a probe of the spec- 
tral properties of a self-adjoint differential operator, say H; it contains "nonlo- 
cal" information not extractible from the much-studied small-time expansion of 
the heat kernel 36, 29 . It reflects the oscillatory flne structure of the eigenvalue 
distribution and, therefore, is directly related to the spectrum of periodic orbits 
[32l [H [I6l [m [20l [HI [m [45] of the nonrelativistic classical-mechanical (or ray- 
optical) system associated with our operator H as quantum Hamiltonian (or wave 
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operator). The interplay among spectral theory, dynamics, and vacuum energy 
is fascinating and rapidly developing, with each subject gaining benefits from the 
others |1UI 151 1431 14UI 1251 135| . Inasmuch as quantum graphs provide instructive 
models of spectral theory and semiclassical dynamics, they should also be commu- 
nicating with vacuum energy. 

Here we show that some elementary techniques recently applied |^ to certain 
traditional boundary-value problems of the Robin type actually apply also to quan- 
tum graphs. Indeed, they may be more valuable there, because of quantum graph 
theory's supply of nontrivial problems susceptible to essentially one-dimensional 
methods. 

Infinite star graphs 

Here we will consider only the simplest type of quantum graph, one consisting 
of a single vertex with N infinite edges attached. (Since many of the issues we will 
study are basically local, there is a sense in which these graphs are the building 
blocks for all others.) The Hilbert space of the model thus consists of vector-valued 
functions, u = {uj{x)} g i^(0, oo)^ , where x as the argument of Uj is the distance 
of the point in question from the vertex along edge j. 




The self-adjoint operator is H = — with certain boundary conditions. We 
impose the usual continuity conditions, 

(1) u,(0) = u(0), Vj = l,...,7V. 
The remaining condition will be one of these: 

• the Dirichlet condition, u(0) = ; 

• the Kirchhoff, or generalized Neumann, condition, 

N 

(2) ^^(0) = 0; 

• our main concern, the Exner-Seba or generalized Robin condition 

N 

(3) ^^(0) = aM(0), q;>0. 

j=i 

Condition Q was apparently introduced in |22| . It is often |21l 138) called the i5 
condition because it can be regarded as the effect of attaching a Dirac delta potential 
at the vertex. At a vertex with only one edge it reduces to the so-called Robin (or 
convective cooling) condition, u'(0) = au{jS). In passing we remark that the label 
"Robin" has almost no historical justification |31j . but it is preferable to "mixed" 
because "mixed boundary conditions" has acquired other meanings |36l I29 | . 



LOCAL DENSITIES NEAR A VERTEX 



3 



Remark. The case a < can also be handled, but the construction of the 
Bondurant transform I^J is then so different as to require a separate discussion. 
With our sign convention, a > is the more "physical" case, where heat flows from 
the hotter to the cooler and H has no negative eigenvalues. 

The Bondurant transform 

In |6] we showed how to obtain solutions of the simplest problems with Robin 
boundary conditions from solutions of the corresponding Dirichlet problems. The 
term "Dirichlet-to-Robin transform" has sometimes been misunderstood as refer- 
ring to an analogue of the Dirichlet-to-Neumann map, hence the justification for 
immediately naming the construction after my junior collaborator. Here one is 
not studying the relation between nonhomogeneous Dirichlet data and nonhomo- 
geneous Neumann data for a fixed solution; instead, one is constructing a new 
solution to a different problem, with homogeneous Robin (or Exner-Seba) data 
replacing Dirichlet data. 

The generalization of the key formula of j6j to an infinite star graph with 
boundary conditions and Q is 

(4) u,{x) ^ iT-^v),{x) 

. N 1 , poo N 

fc=l J •'^ k=l 



1 

a 



Theorem. Ifv{t, x) solves a Dirichlet problem for a constant- coefficient partial 
differential equation [cf. then u{t,x) solves the corresponding Exner-Seba 

problem (though with different initial data). 

Sketch of verification. Q is obtained by solving the ordinary differential 
equation v — Tu, where 

N 

(5) iTu)j (x) = ^ Wfe (x) - auj (x) , 

fc=i 

with the condition of decay as a; — > oo. The heuristics of finding the solution are less 
instructive than the verification that it satisfies all the required conditions. Since 
v{0) = 0, one observes that 

• Uj{0) is independent of j (condition (Q); 

N 

• Tu{x) — v{x), so ^w^(O) — au(0) (condition ||2Jl); 

i=i 

• T and commute with the partial differential operator, so u is still a 
solution. □ 

Integral kernels 

Now consider any of the following initial-value problems. 

(6) Wave: Uu = u^x , u{0, x) = f{x), ut{0, x) = 

(7) Heat: ut = u^x , u{Q,x) = f{x) 

(8) Quantum: iut ~ —u^x = Hu, it(0, x) — f{x) 

(9) Cylinder: uu ^ Hu, m(0, x) = /(a;), u(-foo, 0) = 
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(Of course, there are many other problems involving the operator H on the graph 
that could be considered, but these form a natural and highly useful quartet.) We 
seek the integral kernel (Green function) that solves such a problem via 

(10) Uj{t,x)^Y. Gslit,x,y)fi{y)dy. 

1=1 -^0 

Corollary. Let Gs(t, x,y) be the (matrix) Green function for one of the 
initial-value problems on the graph. Let G{t, \x — y\) be the corresponding 

(scalar) Green function on the real line (also known as the "free" kernel). Then 

(11) Gslit,x,y)^SjG{t,\x~y\) 



Sketch of derivation. The Dirichlet Green function consists of an incident 
term minus an image term, 

GdI (t, X, y) = Sj [G{t, la; - 2/1) - G(t, x + y)] 
— G— — G+ ■ 

In operator language, we want Gs = T^^GdT . (The final T is needed to get the 
correct initial value, SjS{x — y).) In kernel language, therefore, we need 

Gs{t,x,y) =T^^T^GD{t,x,y), 

where f indicates the transpose (real adjoint). But 

(12) (Ttu), 

k 

(13) = ~{Tu)j -2auj . 

From ((T^ and —dyG{t, \x — y\) — +dxG{t, \x — y\) one sees that the incident term 
passes through the similarity transformation unchanged: T~^G-T = G- . From 
(|13|l and —dyG{t, x + y) = -~dxG{t, x + y) one has TyG+ = T^G^ and hence 

(14) Gs}[t,x,y) ^5}[G{t,\x-y\) + G{t, x + y)] + 2aT-][SlG{t, x + y)]. 

(Note that the first two terms of (|14|) solve the true Neumann (not Kirchhoff) 
problem, u'j{0) = Vj.) Working out the last term of (|14|) according to Q), one 
obtains UTl- □ 

The factor (|- - Sj) in ^ wiU come as no surprise to those who are familiar 
with the study of quantum graphs by other methods (see j37j ). 

Remark. The kernel formula Hll|) is correct for a = (the Kirchhoff condi- 
tion), although the intermediate steps are meaningless (in particular, T^^ doesn't 
exist in that case). 

Main results. As corollaries of the corollary, we derive formulas H16(l - (|17|l . 
(j22l, (|23lj I^H)) . PHjl for particular kernels and associated quantities. 
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The wave kernel 

In a one-dimensional system the simplest member of the quartet is the wave 
problem for which the free Green function is (d'Alembert's solution) 

(15) G{t,z) ^ ^[Siz-t)+Siz + t)]. 
Applying and omitting terms that vanish for t > one gets 

(16) Gsht,x,y) - ^Sl[6ix-y~t) + 6ix~y + t)] 



2 

N 



61 ]6ix + y-t) 



a 



where 9 is the unit step function. 

The meaning of (|16l) becomes clearer when one applies (|10|l to get 



(17) u,{t,x) = i[/,(x - t) + f,{x + t)] - y,{t - X) 
1 ^ 



1=1 



N 



e/N 



^ fi{t - X - e) de. 



1=1 



Here we see clearly the incident wave, the immediately reflected and transmitted 
waves from a KirchhofF vertex, and some a-dependent delayed transmission. The 
Robin case (iV = 1) is |5] 



(18) u(t,x) = l[f{x-t) + f{x + t) + f{t~x)]-ae{t-x) 



7(t-x-e)de. 



What is the physical meaning of this time delay? In the context of the wave 
equation, the Robin or Exner-Seba boundary models an ideal spring, or elastic 
support, to which the vibrating medium is attached. The spring absorbs energy 
from the medium and leaks it back out. It is a jolly exercise in integration by parts 
to show that for the solution H18|) the total energy 



(19) 



duV 
'dij 



du 
dx 



dx + -u{t,oy 



is indeed conserved; the field and boundary terms individually are not (unless 
a = 0). 

For later use we note that one integration by parts in 1)19(1 and use of u'{t, 0) — 
au(t, 0) lead to an alternative formula for the total energy, 



(20) 



E : 



du 
'dt 



o u 
dx"^ 



dx, 



in which the boundary term has been formally absorbed into the field term. 



The spectral projection kernel 

Let P{\) be the orthogonal projection operator onto the interval [0, A] of the 
spectral resolution of H . Because H on an infinite star graph has purely absolutely 
continuous (and nonnegative) spectrum, the integral kernel of P(A) may be written 



P{\,x,y) 



cr^LU, X, y) duj, 
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where tr is a well-defined matrix- valued function (not just a distribution in lu). 
Alternatively, a can be defined as the inverse Fourier cosine transform of the wave 
kernel. (It is also the inverse Laplace transform of the heat kernel, and so on for 
all the standard "spectral functions" j36j.) 
That is, we write 

/>oo 

(21) Gs}{t,x,y)= / cos{u}t)asj{uj,x,y)duj 

Jo 

and calculate 



(22) asj{uj,x,y) ^ - coB{ujt)Gsj{t,x,y) dt 



TT 

2 

— 5j sin(cja;) sin(aj?/) 

TT 



2/71 

— „ {Nuj'^ cos[uj{x + y)] + auo sin[cj(a; -I- y)]}. 



In a sense, (|22|l is the ultimate formula concerning the operator H , since all operator 
functions of H can be calculated from it in principle and it contains all facts about 
the spectral resolution in a rather explicit form. 

Kottos and Smilansky |37l Sec. 3B] found the spectral resolution by treating 
the infinite star graph as a scattering problem. For each u) they provide the basis 
of incoming scattering eigenfunctions 



(23) ^/(x) ^ ^/e— 



-<5/ + l(l + e-^'*--^^) 



= -2i5j s\w{ujx) + 2w — ^ e'""^. 

J ' Q.2 _)_ 7V2cj2 

These basis elements are orthonormal (up to the conventional factor ^/2t:) and 
therefore 

1 ^ 

(24) - ^li^Wiy)* = ^5/ (^, X, y). 

1=1 

A calculation verifies that 1)22(1 and 124|l agree. 

Remark: Direct construction of eigenfunctions by applying the Bondurant 
transform to orthonormal eigenfunctions of the Dirichlet problem, while possible, is 
not recommended. In the present problem the immediate results are not orthogonal, 
much less normalized. If a = they are not even linearly independent (because T 
is not invertible), and one basis element needs to be found by a separate argument. 

The local spectral density 

Special interest attaches to the diagonal values of a, 

(25) asj'{i^,x,x) = ^ + ^(-^-1) cos(2Lja;) 



2a /tt 



uj sm(2ujx) ~ — cos{2ujx) 



Clearly here the ^ is the universal Weyl term for a one-dimensional system, the next 
term is the spectral effect of a KirchhofT vertex, and the last term is the Exner-Seba 
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correction. In the limit a — > +00 there is some cancellation between the second 
and third terms, resulting in 



(26) 



cos(2a;a;) — (T_d/ (w, x, x). 

TT TT 



as expected for a Dirichlet vertex. 

Because the spectrum is continuous, one can't integrate to get a density 
of states. However, subtracting off the Weyl term and paying due attention to 
distributional limits, one can obtain a meaningful global spectral density: 

N r 



(27) 



Ap{u;) 



E 

N 



dx 



5{u) 



Na/n 
2 + Ar2^2 



■5{lo) 



This expression approximates the incremental effect that such a vertex would have 
in a problem with discrete spectrum. The first term in H27|l is the Kirchhoff term; it 
vanishes when N = 2, because a Kirchhoff vertex with exactly two edges is vacuous. 
The other term is the Exner-Seba correction; it vanishes when a = because its 
first term distributionally approaches ^S{uj) in that limit. Alternatively, 127|l can 
be simplified to 



(28) 



N 

Ap{^) ^ - - Siiu) 



a' 



Na/n 



here the first term is the correct formula for a Dirichlet vertex and the remaining 
term is 0{a~^) as a — > 00. 

Remark: The meaning of a Dirac delta distribution in formulas such as H27(l 
and 1)28(1 is that the spectral counting function N{ijj) has a nonzero limit as uj 
approaches from above {N{uj) being understood to be for negative uj). For 
example, (I28|l is simply the derivative of the formula 



(29) 



AN{u) 



N 

T 



tan 



Nuj 



for the incremental effect of the vertex on the total number of eigenvalues in the 
interval < A < w^. 

The Bondurant method cannot be applied directly to a finite interval, be- 
cause no transformation T will work for both boundary conditions simultaneously. 
However, the appropriate operators for the two boundaries can be applied al- 
ternately to construct a solution as an infinite series (a sum over closed classical 
paths, generalizing the classic method of images). In [S] the wave kernel and hence 
the local and global spectral densities were obtained in that way for a finite inter- 
val with one Robin and one Dirichlet endpoint. Numerical evaluation reveals the 
correct eigenvalues for the problem emerging as spikes in the global density (the 
counterpart of ((23 )■ See |44) for a related study in two dimensions. It should be 
straightforward to extend this analysis (and also the study of vacuum energy, etc.) 
to an arbitrary finite star graph, and in principle to more complicated quantum 
graphs. It is noteworthy that in these systems no semiclassical (or stationary-phase) 
approximation is needed to obtain the representation of the spectrum in terms of 
classical paths; the only approximation involved is the truncation of the sum at 
some maximum path length if and when one resorts to numerics. 
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Heat and quantum kernels 

The same machine lfTT|l can be used to treat the problems and for which 
the free kernels are 

G{t, z) = (4^t)-i/2e-.V4t^ z) = {Amt)-^l^e~^"/^'\ 

respectively. (The results are qualitatively similar to (|35|l below, with the com- 
plementary error function appearing instead of the exponential integral function.) 
Studying the heat kernel is the traditional route to equations like H29|l for par- 
tial differential operators. In one-dimensional systems such as quantum graphs, 
however, the wave kernel built from p5(l appears to be easier to calculate. 

Using the heat kernel, the Robin case of H28|l was obtained in |S| Sees. 3.3 and 
5.5]. That analysis extends to flat boundaries with constant a in any dimension. 
The results fit naturally with those of [J in dimension 3 and |44| in dimension 2, 
where the leading orders in boundary curvature are also included. All these formulas 
are exact in a; of course, when the heat-kernel formulas are expanded in power 
series in a they match and extend the relevant terms tabulated in such references 
as 36, 29j. (See also [191 |7]-) Apart from a unifying point of view, it is not 
claimed that these results are particularly new; in fact, the Robin heat kernels in 
dimensions 1 and 3 were found in 1891 by a closely related method jllL I12| . 



Vacuum energy density 

First, a paragraph which pure mathematicians are free to ignore: From a 
physical point of view, vacuum energy involves a relativistic (usually massless and 
bosonic) field. (It is of no relevance, therefore, to quantum graphs if they are re- 
garded solely as models of nonrelativistic electrons in networks of wires.) Formally, 
the total energy corresponding to the wave operator H is 

(30) ^ " 2 ^ ^" " 2 / ^ ^^^^ 

n=l •'^ 

for an operator with purely point spectrum. (To avoid irrelevant complications, let 
us also always assume that H has no negative spectrum.) Equally formally, the 
local energy density is 

1 

(31) Too{x) ^ - / uj(T{uj,x,x)dLj 

^ Jo 

(without the requirement of point spectrum). The origin of these expressions is 
the "second-quantized" theory of a field satisfying the wave equation © , in which 
each normal mode of the field (with frequency lu) becomes a quantum harmonic 
oscillator (with ground state energy iw). Then (|31|l results from the integrand of 
()20|l . and 1)3011 comes from integrating (|31l) over all space or just from adding up 
the energies of all the modes. Both integrals, (|3())l and (|31)l . are divergent at the 
upper limit and are to be defined by a renormalization procedure. 

Our claim is that vacuum energy should be of mathematical interest even in 
models whose physical relevance is questionable. Therefore, we provide here a 
precise mathematical definition, which incorporates a particular renormalization 
prescription (whose physical rationale need not concern us): Consider the cylinder 
kernel (the Green function of ^) on diagonal (y = x, I — j), find its asymp- 
totic expansion as t ^ 0, and extract the coefficient of the term proportional to t, 
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times — ^; this is the vacuum energy density, Toq{x). When appropriate, integrate 
over X (and sum over j in our graph problem), before taking t to 0, to define a total 
energy, E. 

The intuition behind this definition is the following. Let G{t, x, y) be the cylin- 
der kernel of the problem under study. (For an infinite star graph it is the ma- 
trix Gs .) Then 



(32) G{t,x,x)= e-'^'a{LU,x,x)duj 







and (when "appropriate" ) its trace is 

r°° . °° 

(33) TTG{t)= ^G/(i,x,x)da; = ^e"""*. 

•^0 j=i ji=i 

Now take — ^ ^ of (|32|l and H33|) and let t approach 0, formally obtaining H31|) and 
(|30|l . respectively. Systematically throwing away the terms of negative order in the 
Laurent expansions (and a possible logarithmic term), one arrives at our definition. 

To find the vacuum energy of an infinite star graph, one can apply the Bon- 
durant machine one more time, to the free cylinder kernel 

(34) G{t,z)^-^, 



obtaining 



<2 + (x - y] 

2a 



^^^^ ~ t^ + ix- yf + U " Jt^ + {x + yf 



+ ^e"(^+^)/^Im 



»*/^ Ei - ^ (X + y)") 



It follows that 

The most important parameter in this problem is the dimensionless product 
ax. Therefore, at short distance the Kirchhoff term dominates: 

(37) + + ^ In \ax\ + 0{a^x\ 
whereas at large distance the energy density is almost pure Dirichlet: 

(38) Tooix)^ -^ + 0{a-'x~''). 

The nonintegrable 0{x^^) singularity in (|37l) would interfere with calculating a 
total energy by integration over x, even if the edges were finite. The renormalization 
procedure implicit in our definition does not commute with the spatial integration, 
however, and it leads to a finite total energy [SI 1241 As agreed, we will not delve 
here into the physical issues thereby raised (which are still somewhat controversial) . 

In a sense the calculation based on the cylinder kernel was unnecessary, given 
the spectral formulas (|22|l . (|25() . and H27() . Indeed, H36() can be obtained from 
(I32|l and 1)2 5|) and a Laurent expansion, or even immediately from (I31II using the 
subtracted spectral density appearing in the integrand of l|77l) . (The subtraction of 
the Weyl term corresponds to the subtraction of the leading Laurent term. In other 
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problems, such as higher dimensions, additional subtractions would be necessary.) 
Similar methods were used in |42| to find the effect of a flat Robin boundary in 
any dimension. In general, however, it is easier to find the small-i expansion of a 
cylinder kernel than to find the detailed spectral resolution; as for heat kernels, one 
expects useful calculations to run in the opposite direction. Romeo and Saharian 
|42j also gave complicated integral formulas for the vacuum energy and energy 
density of a finite interval with two Robin boundaries. Our methods will instead 
give these quantities as infinite sums over classical paths |27| . 

A broader perspective 

Now let if be a generic differential operator (self-adjoint, elliptic, positive, 
second-order, with scalar principal symbol) in dimension d. Let T be the cylinder 
kernel and K be the heat kernel. Then 



TtT= e-^'^dN^ e-*'^p{u;)duj, T{t,x,x)= e-*'^cr(w, x, x) dw, 
Jq Jq Jo 

/•OC /"OO /"OO 

TrK^ e-*^dN= e~*'^%{Lu) du , K{t,x,x)= e-*'^\{uj,x,x) duj, 



where N now is the number of eigenvalues less than or equal to A = o;^. 

It is known that as i — > the global quantities have expansions of the forms 

oo oo 

(39) Tr T - ^ e,t-^+^ + ^ In t, 

s~d odd 

oo 

(40) TtK ^J2b,t^-'^+'^/\ 

The local quantities, T{t,x,x) and K{t,x,x), have (nonuniform in x) expansions 
of precisely the same respective forms; we do not introduce a separate notation for 
their coefficients. If d — s is even or positive, 

(41) e, = 7r-i/22'^"*r((d-s + l)/2)6,. 
If d — s is odd and negative, then 

^_l^is-d+l)/22d-s+l 

^^'^^ V^T{{s-d+l)/2) ^ 

but, most strikingly, in that case 

(43) Cs is undetermined by the br 
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I have expounded in detail elsewhere "26", "24", "25' how the connections H41|l - 
(j43|l come about: The bg are proportional to coefficients in the high-frequency 
asymptotics of the Riesz means 33j I34| of A'^ (or of / a) with respect to A, while 
the Cs and fs are proportional to coefficients in the asymptotics of the Riesz means 
with respect to oj. Alternatively, they are related through the poles of the zeta 
functions of g and [28ll30l l4]. 

The new half of the cylinder- kernel coefficients (those in — of which the 
first, ed+i , is proportional to the vacuum energy — are a new set of moments of 
the spectral distribution. What are they good for, mathematically? Unlike the old 
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ones, they are nonlocal in their dependence on the geometry of the domain and the 
coefficient functions of H. They probe, in a comparatively crude, averaged way, 
the oscihatory spectral structures that are correlated more precisely with periodic 
and closed classical paths. Thus they embody, at least partially, the global dy- 
namical structure of the system; they are a half-way house between the heat-kernel 
coefficients and a full semiclassical closed-orbit analysis. 
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